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ABSTRACT 

We discuss short wavelength (inertial wave) instabihties present in the standard two- 
fluid neutron star model when there is sufficient relative flow along the superfluid 
neutron vortex array. We demonstrate that these instabilities may be triggered in 
precessing neutron stars, since the angular velocity vectors of the neutron and pro- 
ton fluids are misaligned during precession. The presence of such an instability would 
render the standard, solid body rotation, model for free precession inconsistent. Our 
results suggest that the standard (Eulcrian) slow precession that results for weak drag 
between the vortices and the charged fluid (protons and electrons) is not seriously 
constrained by the existence of the instability. In contrast, the fast precession, which 
results when vortices are strongly coupled to the charged component, is generally un- 
stable. This implies that fast precession may not be realised in astrophysical systems. 



1 INTRODUCTION 

Over the last few years, the modelling of neutron star 
free precession has attracted considerable attention from 
' theorists. This has been caused in large part by the 
\ report ed observations of free prece ssion in PSR B1828- 
. 11 by IStairs. Lvne fc Shemaj lj2000l ). There exist several 
' other precession candidates, mos t notably PSR B 1642-03 
. (|Shabanova. Lvne fc Uramall200ih . but none are as convinc- 
■ ing as PSR B1828-11, whose pulse variations extend over 
' multiple precession periods and are visible in multiple chan- 
, nels (pulse time-of-arrival and beam shape). 

This precessional motion is of great interest in several 
, regards. Firstly, PSR B1828-11 is the only clear precession 
candidate out of the several thousand known pulsars, which 
begs the questions why is precession so rare, and what has 
set PSR B1828-11 into precession in the first place? Sec- 
ondly, the free precession timescale as extracted from the 
observations is approximately 500 days (modulo a possi- 
ble factor of two up or down, depending on which of the 
observed harmonics in the frequency spectrum is the funda- 
mental one) . This is extremely uncomfortable from the theo- 
retical point of view, as the canonical model of a neutron star 
containing one or more superfluid components would suggest 
a much shorter precession timescale. Instead, this long pe- 
riod free precession period seems to be consistent with the 
simple wobble of a star deformed away from sphericity by 
strains in its crust or by a ma gnetic field, with superfluid- 
ity p laying no role whatsoever l|Akgun. Link fc WassermanI 
I2OO6I ). It is this second question regarding the apparent lack 
of the effects of superfluidity with which we will mainly be 
concerned in this paper. 

The basic picture is as follows. Both on theoretical 
grounds, and in order to explain the glitch phenomena, neu- 
tron stars are believed to contain superfluid components 



whose rotation is achieved by their forming arrays of vor- 
tices poi n ting a long the rotation axis. In a pioneering paper 
IShahanJ (|l977t ) showed that in the fimit where these vor- 
tices attach or 'pin' themselves tightly to the radio-emitting 
crustal phase (as is required to explain the glitches) , the clas- 
sical free precession of the star is greatly modified by the gy- 
roscopic effect of the superfluid, whose angular momentum 
vector would be fixed in the crust's frame. The resulting 
free precession period is then approximately P/{Is/Ii + e), 
where P is the spin period. Is and Ji the moments of inertia 
of the superfluid and crustal phases, and e the dimensionless 
non-sphericity in the crust's moment of inertia tensor. Note 
that by 'crust' we really mean the crystalline crust itself and 
any other part of the star coupled to it on sufficiently short 
timescale to participate in its precession, which probably 
includes the core proton fluid. 

One region where neutron vortices are expected to in- 
teract strongly with the crust is the inner crust, where neu- 
tron and crustal matter coexist. The interaction here is ex- 
pected to be sufficiently strong that the vortices pin tightly 
to the crustal lattice, resulting in a very short free pre- 
cession period, in conflict with the PSR 1828-11 observa- 
tio ns. A possible r esolut ion to this problem was suggested 
by iLink fc Cutleil (|2002l l. who argued that forces created 
by the precessional motion itself would un-pin the vortices, 
restoring long period free precession, consistent with the 
observations. More problematically, neutron vortices in the 
core of the neutron star are also expected to interact with 
the charged component. SpecificaUy, the interior magnetic 
field is composed of a large number of fluxtubes (a conse- 
quence of the protons being condensed in a type-II super- 
conductor), which can interact with the vortices via a pro- 
cess known as mutual friction. As was shown in detail by 
ISedrakian. Wasserman fc Corded l|l999l ), even if this con- 
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pling force isn't strong enough to completely pin the vor- 
tices, the free precession timescale in the coupled system 
is likely to be much shorter than the observed precession 
timescale. The s ignific ance of this result was emphasised re- 
cently by iLinkI (|2003l ). who argued that the observed slow 
precession necessitated a profound change in our view of 
neutron star interiors, requiring either that one or other (or 
both) of the posited superfluid/superconducting phases do 
not exist, or else that there are no regions in the interior 
where they exist simultaneously, perhaps because the pro- 
tons form a type-I superconductor. 

It is this line of argument we address in this paper. 
The results outlined above, leading to the prediction of fast 
precession, have all relied implicitly on treating the super- 
fluid and non-superfluid phases as rigid bodies, each with 
a well defined angular velocity vector. However, in reality 
the phases are better described by the fluid Euler equa- 
tions, the ansatz of rigidity having been made mainly for 
reasons of tractability. In this paper, which expands on ini- 
tial re sults discussed by iGlampedakis. Andersson fc JonesI 
l|200i1 ). we demonstrate that the relative motion between 
the two phases that rigid body free precession generates can 
give rise to a superfluid two-stream instability. The instabil- 
ity itself occurs in the inertial modes of the rotating fluids, 
and is related to the experimentally observed instability of 
Heli um II explained long ago by Glaberson and collabora- 
tors (|Glab erson. Johnson fc Ostermeieil[l974l ). It is also in- 
timately lin ked to the neutron star instab i hty d escribed very 
recently by ISiderv. Andersson fc Conierl (120071 ') . As we will 
argue in detail, this instability renders the free precession 
rigid rotation ansatz inconsistent and makes all subsequent 
conclusions concerning the interior composition of neutron 
stars unsafe. 

The plan of this paper is as follows. In Section 2 we de- 
scribe the basic formulation of the two-fluid coupled system. 
Section 3 provides the plane-wave analysis and the instabil- 
ity criteria for inertial modes in various regimes. In Section 
4 we discuss the relevance of these results for freely precess- 
ing neutron stars. The detailed precession modes are anal- 
ys ed in Appendix A drawing heavily on t he previous work 
of lSedrakian. Wasserman fc Corded l| 19991 ). Our conclusions 
can be found in Section 5. 



2 FORMULATION 

When viewed as a multifluid system, a neutron star core 
can be modelled (at the simplest level) as a mixture of su- 
perfluid neutrons and a conglomerate of superconducting 
protons and normal electrons. In this picture it is assumed 
that the charged components are electromagnetically cou- 
pled on a timescale which is short compared to the dynamics 
that is being investigated. The charged plasma then forms 
a charge neutral fluid, which we will loosely refer to as the 
"protons" in the following. When required we will use the 
constituent index x, which can be either n or p, to iden- 
tify the different components. A key property of the sys- 
tem is that the neutron and "proton" fluids do not flow 
independently with respect to each other. One of the key 
coupling mechanisms is the vortex mediated mutual fric- 
tion, which acts dissipatively to prevent relative flow and 
also affects the nature of the waves in the system. The 



mutual friction force appears as the result of interactions 
between the neutron vortices and the charged component. 
In the standard form, it represents scattering of electrons 
off of the magnetic fie l d ass o ciated with the vortex core 
jAlpar. Langer fc Saulsl Il984l : lAndersson. Siderv fc Cornell 
l2006h . but there may also be contributions due to the in- 
teraction between vortices and magnetic fluxtubes in the 
charged fluid etcetera. This latter effect may be important 
since the proton s are expected to condense i nto a type II 
superconductor (|Bavm. Pethick fc Pines|[l969l ') in which the 
magnetic flux is carried by fluxtubes. 

The s mooth-averaged hydrodynamical equ ations for our 
system are ij Andersson. Siderv fc Comeijl2006l ). in a rotating 
frame with angular velocity S7o, 

{dt + viV,)v7 = V.-^n + 2e,,kviQ.'^, + ff (1) 

{dt + vl,V,)v^^ = V.Vp + 2e,,fcW^nS - — /r' + i/oo V'«; (2) 

In these equations are the constituent velocities and 
Xp — Pp/pn is the density ratio. Since it is expected that 
the proton to neutron ratio in the outer neutron star core 
where the above equations apply is about 0.05, Xp may also 
be taken to represent the proton fraction Pp/(pn -I- Pp). We 
will not be making a distinction between these quantities in 
the following. In the two Euler equations, the scalar func- 
tions Vx = Mx + "J" represent the sums of speciflc chemical 
potentials and the gravitational potential. Since our focus 
will be on inertial waves, we can assume that the two fluids 
are incompressible. This means that we also have the two 
continuity equations 

V.wi = (3) 

In writing the above Euler equations we have made a 
number of simplifying assumptions. We have ignored the en- 
trainment effect between neutrons and proton^ . We have 
assumed that the charged components are comoving, which 
when combined with an assumption of local charge neutral- 
ity, eliminates the electromagnetic Lorentz force from the 
equations of motion. We are also not considering the effect 
of vortex tension or the elasticity of the vortex array, even 
though quantised neutron vortices are present (in the vol- 
ume averaged sense) in ([1]). The effect of the tension on the 
kinds of waves that we will study has been considered by 
I Siderv. Andersson fc Comer! l|2007h . who found that it was 
generally safe to ignore it. We are also not accounting for 
the analogous fluxtube tension in ((2)1 . 

We are, however, including the effect of shear viscos- 
ity in the equations; this is important since we will con- 
sider very short wavelength dynamics. In a superfluid neu- 
tron star core, electron-electron scattering pr ovides the dom- 
inant contribution to the shear viscositv dFlpwers fc Itohl 
1 19761 : lAndersson. Comer fc Glampedakis! l2005l ). This leads 
to the last term on the right-hand side of It should 

^ Even though the entralnment, due to the strong interaction be- 
tween neutrons and protons, is an important effect we ignore it 
in order to simplify the analysis. Its inclusion would affect the 
quantitative results, but should not lead to any qualitative dif- 
ferences. It would, of course, be interesting to study the role of 
entrainment in this problem at some stage. 
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be noted that this description is expected to be accurate 
at temperatures significantly below the various superfiuid 
transition temperatures. At higher temperatures one would 
also have to account for various excitations, which will lead 
to the presence of a numbe r of additional viscous terms (see 
lAndersson fc Comeij (j2006h for a discussion) . In order to es- 
timate th e kinematic viscosity coefficient in ^ we use the 
results of lAndersson. Comer fc Glampedakid ( 2005f ). For a 
uniform density star with mass 1.4Mq and radius 10 km 
(our canonical values) we then find that 



10^ 



l(f_K 

T 



2 / 

cm /s 



(4) 



The mutual friction force in ^ and dH) has the stan- 
dard form (|Hall fc Vinen|[l956l : lAndersson. Siderv fc Conierl 
l2006l ) 



/■mf f2 

Ji — oeijk€ 



kml ^7 n np 



np I yj! 1 k 



(5) 



where = 2f2Jj -|- e*-* Vjv'j^ is the neutron vorticity, the 
relative fiow is given by Wnp — v^ — Vp and a 'hat' denotes a 
unit vector. This form for the mutual friction follows from 
assuming a resistive drag force (per unit length) 



(6) 



acting on each neutron vortex (moving with velocity vl^). 
Balancing the drag force against the standard Magnus force 
acting on the vortex leads to ([5|. The dimensionless coupling 
coefficients B and B' then follow from 



B 



n 



l + 7^2' 



and 



B' = 



1 + 7^2 



(7) 



where TZ — TZ/ pjiH. Here k = /i/2mn is the quantum of cir- 
culation. This description is quite generic. It can be taken to 
account for dissipative scattering of electrons off the mag- 
netic field of the vortex core . As d iscussed by, for example, 
lAndersson. Siderv fc Comeij l|2006D this leads to a typical 
value of B « 4 X 10~*. This means that we have 



B ^ 7^ < 1, S' ^ B^ < B 



(8) 



We will refer to this as the weak drag limit. 

However, it may well be the opposite limit that ap- 
plies. Strong coupling between the neutron and proton fiu- 
ids may originate from the interaction of v ortices with 
the much more numerous magnetic fluxtubes l|Saulslll989l : 
iRuderman. Zhu fc Cheiilll998l : iLinkI 120031 ). In this picture, 
each vortex is "pinned" to a large number of fluxtubes. As 
a consequence it tends to be forced to move with the pro- 
ton fluid. Essentially, this additional interaction can be ac- 
counted for by modifying the drag force to 



/d = 7^(Wp ~ vl) + n'{vl - vl) 



(9) 



where i)^ is the velocity of the fluxtubes. If we assume that 
the fluxtubes move with the electron/proton fluid, then it 
is clear that we will again end up with a mutual friction 
force of the form ((S)). Of course, this will not be the case if 
there is a slippage between the fluxtubes and the charged 
componen10. We will not consider this possibility here. The 
precise form of the drag force experienced by the vortices in 



^ Neither will ((5J result if the drag coefficient TZ is velocity de- 
pendent. 



this case is not well known, but it is usually assumed that 
the pinning to the fluxtubes is efficient. This would mean 
that one would have TZ' ^ TZ and the final mutual friction 
force takes the form ([Sjl with TZ ~ TZ' /pnK S> 1. In this case 
we have the strong drag limit, 

B^^<1, 1-B'^B''<1 (10) 

It is worth noting that some constraints on the value 
of B can be placed by post-glitch relaxation data, using 
the two-fiuid formalism (with the crust assumed rigidly at- 
tached to the "protons"). For example, Vela's 1988 "Christ- 
mas" glitch i mplies a crust-core coupling tim escale Tc < 
300(7i//s) s l|Abnev. Epstein fc Qlint"ol 1 19961 ). A simple 
glitch relaxation model can be constructed from eqns. ([l}- 

assuming rigid-body rotation for the fiuids. The result- 
ing coupling timescale due to mutual friction is found to 
be Tmf = Is)/(2Q.B), in agree ment with the result of 
lAndersson. Siderv fc Cornell (|2006l ). For a rotation period 
P = 27r/n~0.1s and 7i « 0.17s, the requirement Tmf < Tc 
provides the constraint B > 2 x 10~^. This translates to 
either 7^ > 2 x 10"^ (weak drag) or 7^ < 4 x lO" (strong 
drag). 

Finally, it is worth making a few remarks on 
the precession model that we co nsider. Following 
ISedrakian. Wasserman fc Corded (|l999l ) we assume that 
the global motion follows from ([T]) and ^ after impos- 
ing that each fiuid rotates as a solid body. Misaligning 
the two rotation vectors we then arrive at a convenient 
two-component model for free precession which accounts 
for the damping due to mutual friction. The detailed 
precession modes are discussed in Appendix A. This model 
represents the current state-of-the-art in this problem area. 
However, it is obviously simplistic. The true "layering" 
of a real neutron star, and possible effects due to the 
various interfaces/boundaries are not accounted for. This 
means that one must be careful when trying to connect 
our model to detailed calculations for the various param- 
eters, e.g. the mutual fric ti on co efficients discussed by 
lAndersson. Siderv fc Cornell l|2006l ). Basically, the model 
should be seen as a body averaged representation of the real 
system, and the same is true for (say) the drag coefficient 
TZ. 



3 LOCAL WAVE ANALYSIS 

The present investigation is motivated by the fact 
that a relative flow along a superfluid vortex ar- 
ray may trigger an instability that leads to the 
formation of vortex tangles and a state of super - 
fiuid turbulence, see lAndersson. Siderv fc Comer! l|2007h : 
ISiderv. Andersson fc Comerl ( 2007^ for discussion and ref- 
erences to the relevant literature. Our aim is to find out 
whether this vortex instability may be active in a precessing 
neutron star. 

As a first step we consider perturbations with respect 
to a stationary background configuration where both fluids 
rotate rigidly, i.e 

e^-'^{Q^ — Q,Oj)xk, — constant (11) 



allowing for an arbitrary orientation of the angular velocity 
vectors. After linearising and flxing f^o = fi^ = finW*, i.e. 
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working in the neutron frame (for later convenience), the 
various quantities appearing in the Euler equations ©-((S]) 
become 



= V'xO + 5ip^ 

We now focus our attention on short wavelength motion, 
such that we can analyse the dynamics without consider- 
ing boundary conditions etcetera. In practice, we use the 
standard plane-wave decomposition 



constant 



(12) 



in IT]), ([2]) and ([3]). It is useful to decompose various vectors 
into a piece parallel to the vortex array and a piece that is 
orthogonal to it, i.e. use 



— A^ h -f A^ 



(13) 



and similarly for all other variables. Given that we have 
assumed that the fluids are incompressible, we must have 



(14) 



In other words, the waves we consider are transverse. 
To simplify the analysis somewhat, we will assume that 
the waves propagate along the vortex array in the 
following. This means that we align fc' with n\ see 
ISiderv. Andersson fc Cornell (|2007l ) for comments on this as- 
sumption. Hence we have 



fc' = kufi 



and 



hjAi 







Al 







(15) 



(16) 



Finally, we define two additional unit vectors and A* to 
form an orthonormal triad. 



s^fii = 



and 



A = e SjUk 



(17) 



One of the key reasons why the above assumptions are 
useful is that we do not need to worry about the pertur- 
bations in the potentials i/>x. They are simply obtained by 
projecting the Euler equations onto n*. To determine a dis- 
persion relation for the propagation and analyse the nature 
of the associated waves, it is sufficient to consider the equa- 
tions obtained by projecting the Euler equations in the di- 
rections and A*. 

Finally, in order to write down the linearised version of 
the Euler equations ([1} and ([2} we need an expression for 
the perturbed mutual friction force. This can be written 



Sfi^t = -2QnB{Al - Al) 
+ ik\ 

+ 2^,-,B'^'''hj{vl - vl) + ifc||B'™o,(n-'Aj, - fCAi) (18) 



+ ifc||B ejkiW'Qn^ A^^ + m||e''''°nj Afe^ 



3.1 Restricted problem: perturbations in neutron 
fluid only 

Before considering the general solutions to the above 
equations it is useful to note that a simplified ver- 
sion of the present problem has alre ad y bee n dis- 
cussed bv ISiderv. Andersson fc Comer! (|2007l ') (see 



IClaberson. Johnson fc Ostermeieil (1974) for the anal- 
ysis of the superfluid Helium case). They assume that 
the proton fluid is locked to its background value, 



i}'j)xk in our case, and that it is sufficient 



to consider only the neutron equation ([l|. It is useful to 
recall the results obtained for this simplified problem for 
two reasons. First of all, the discussion explains our general 
strategy. Secondly, and more importantly, we will later 
demonstrate that the results remain essentially unchanged 
in the weak drag limit. 

Assuming that the protons remain unperturbed, we eas- 
ily obtain 

{ia + 2nnB - iB'kiiwii) Ai 

+ [ (1 - B')2nn - iSfciiwii ] e^'^hjAl = (19) 

where 

W\\ = n^{VnO - Vpo) = -"-iWpO (20) 

Taking the inner product of p9|) with v^, gives 

{Alf = ^ s,A:, = ±i{kAl) (21) 

Hence the modes under discussion represent circularly po- 
larised waves. Next, projecting (|19|l along s* and A' we get 



[2Bnn+j(cr-S'fc,|™||)] {s.Ai,) 

+ [2{1 ~ B')Qn - iBw^i fc|| ] ikAl) = (22) 

[2(1 -B')nn-iBwfy] (S.Ai,) 

- [ 2Bnn B'k^^wn )] {X,Al) = (23) 

Requiring a non-trivial solution for this system we find the 
dispersion relation, 

a= [S'fc||m|| ±2(1-S')!^n] + (2!^n =F ) (24) 

Note that this result is valid for arbitrary B and B' , i.e. 
it is relevant for both the strong and weak drag lim- 
its. It is also in agreement with the result found by 
IClaberson. Johnson fc Ostermeierl lll974|) in the limit of 



purely transverse waves, and Siderv. Andersson fc Comer! 
(2007) in the limit B' 0. 

In the limit of vanishing mutual friction the dispersion 
relation reduces to cr = ±2f2n. These modes represent col- 
lective vortex Kelvin oscillations. At the macroscopic level 
these wave s can b e identified with incr tial waves in the neu - 
tron fiuid (jSonin! [1987 : Siderv. Andersson fc Comei! [2OO7I ') . 
In the presence of mutual friction the inertial modes are 
generally damped but, remarkably, they may become unsta- 
ble when there is a relative flow along the vortex array. We 
have an instability when 

20^ 



Im CT < 



hill > 



(25) 



In the case of superfluid Helium this instabili t y was flrst 
discussed bv !Glaberson. Johnson fc Ostermeiei! 1 1974h . and 
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is o ften refereed to as the Donnelly- Glaberson instabil- 
ity (jBarenghi. Donnelly fc VinenI l200ll ). In the context of 
neutron stars the instabili t y was recently discussed by 
ISidery. Andersson fc Cornell l|2007h . who argued that is a 
typical "two-strea m" instability, a generic property of multi- 
fluid systems (see lAndersson. Comer fc Prij ()2004l ) for ref- 
erences to the relevant literature). The relevance of this 
instability fo r neutron s tars has also been discussed by 
IPeralta et all (|2005l . l2006l ). 

We want to consider these plane- wave results in the con- 
text of a neutron star that is freely precessing. This means 
that we need to connect the background that we have used in 
the plane- wave analysis to the global motion of a precessing 
star. In all current models of precession (see the discussion 
in Appendix |AJ it is assumed that the two components ro- 
tate as solid bodies. However, since the rotation axes are not 
aligned one would in general expect both fJ^ and Q,p to be 
time dependent. Since they are taken to be constant in the 
plane-wave calculation it is obvious that the result is only 
valid on a timescale that is short compared to the variation 
in and fip. As we will argue later, the relevant timescale 
to compare to is the precession period. Hence our analysis is 
applicable only if the instability growth time is significantly 
shorter than the precession period. 

It is also clear that the unstable waves need to have rela- 
tively short wavelength. This becomes obvious if we consider 
the explicit form for wu. For our problem we have 



[mmx^) - [m){>^^^') (26) 



This shows that the projection, w\\, of the background flow 
along the direction of the neutron vortices varies with po- 
sition. This obviously does not invalidate the plane-wave 
analysis, but it is clear that the characteristic lengthscale 
must be much smaller than the radius of the star. 

Finally, we need to appreciate that there will 
also be a cut-off at short wavelengths. If we had 
accounted for the vortex tension in our analysis 
we would hav e found a slightly different i n stabil- 
ity criterion (Glaberson. Johnson fc Ostermeieij 1 19741 : 



ISidery. Andersson fc Comeij|2007l ) 



2nn 
1^ 



+ vk\\ 



(27) 



where u represents the tension. From this it is clear that the 
tension tends to stabilise the system for large k\\ (short wave- 
lengths). However, if we compare the relative magnitude of 
the terms in (I27|l we immediately see that the tension is 
rather insignificant. The two contributions are equal for a 
wavelength. 



Ao = 27rf^J ^ 7 X IQ-^ s/Pis) cm (28) 

where P{s) is the rotation period in seconds. In other words, 
the tension contribution dominates only for wavelengths Ao 
comparable to the intervortex spacing 



~ 10"V'P(s) cm (29) 
At that level the macroscopic hydrodynamics that we have 



been using is no longer relevant. Equations (O-® arrive af- 
ter averaging over lengthscales much larger than dn. Conse- 
quently, the plane wave analysis makes sense provided that 
A S> dn, which also means that A 3> Aq. Hence, as far as 
the present hydrodynamics problem is concerned, we can 
always neglect the vortex tension. The vortex tension would 
be important if we were to consider the dynamics of a sin- 
gle vortex. In fact, an analysis of that problem shows that 
an oscillating vortex can also become unstable. The relevant 
instability criterion is similar to l|27|) . although as expected 
for large fcy the first term is irrelevant and the tension plays 
the key role. 



3.2 Complete problem: perturbations in both 
fluids 

In a neutron star core there is no compelling argument for 
clamping the proton fluid as we did in the previous example. 
It was a useful assumption because it simplified the algebra. 
It also provided a straightforward demonstration of the pres- 
ence of unstable waves in the system. However, if we want to 
make contact with realistic precessing neutron star models 
we need a more detailed analysis. Hence, we consider the 
case where both fluids are perturbed. 

In the general case we have the two equations of motion 



2n„BA^ + 2f7nS'e'^*Ti,^P 



and 



icr-|-2rinS-iS'if||fc||)X 
- 2n„BA^ + 2 
+ [20.^(1 - B') - iSwii/fcii] e'^fijAl = (30) 



+ 



ik\\W\i 



2nn 



+ 



ik\\w\\ 



B 



+ 



P 

t^n ( 1 



+ 20. 
2B 



B' 



^^%Al^O (31) 



where Qp — hiQp. Projection of these equations along the 
and A' vectors results in the 4x4 system 

{ia + 2Bfln - iB'wfy) {siAl) - 20nB{siAl) 

+ [20n(l - B') - iBwfy] (X^Al) + 2!:!„S'(AMp) = 

(32) 

{ia + 2Bnn - iB'«;||fc||) {\,Ai) - 2nnB{X^Ai) 

- [2Qn{l - B') - iSwiifcii] (s>A^) - 2nnB'(s,4) = 

(33) 



ia + Vcckl + 2riii ( — 



ik\\W\\ 



20^ 



B' 



- ik\\w\\ 



20^ 



B' 



^T7 



+ ik\\W\\ 



B_ 



(5.AL) 

(A,A^) 

(A,yi;)=0 (34) 
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ia + Ucckfi + 2Sln ( — 



- ik\\W\\ 



+ 



■I. I ^' 

ifc||ti)|i — 



20, — 



2nn ( — 



B' 

1 - 2 — 



B 



{s,Al)=0 (35) 



The dispersion relation associated with this system is a 
fourth order polynomial in a. The roots can be explicitly cal- 
culated, but the dependence on the various parameters will 
be complicated. In order to explore the nature of the waves 
in the system it makes sense to focus on specific limiting 
cases. This way we also obtain useful approximate expres- 
sions for the growth time of the unstable modes. As a first 
approximation it is natural to consider the inviscid problem, 
i.e. we set v^c = 0. 



3.2.1 Approximate inviscid solutions 

With the shear viscosity eliminated, we can arrive at differ- 
ent approximate solutions by considering, separately, (i) the 
short wavelength limit fcyii — > oo and (ii) the weak/strong 
drag hmits, 7?. — > and TZ —> oo, respectively. 

The short wavelength approximation leads to the fol- 
lowing mode solutions, valid for every TZ, 

0-1,2 ~ ±2nn + [iB =F B') {2fln T (36) 

II 20 
asA ~ ±{nn + ^1) + + —{iB T B') (37) 

Xp 

The top pair is identical to the modes found in the problem 
where the protons were kept clamped. These modes can still 
become unstable under the action of mutual friction, and the 
instability criterion (|25p remains unchanged. 

For a generic value of TZ both mode pairs in (|36p - (|37p 
describe waves where the neutron and proton fluids oscillate 
in concert. However, a closer look reveals their distinct na- 
ture. In the weak drag limit, TZ ^ 0, equations pO|l and (|3ip 
decouple, and the modes (I36|) represent pure neutron fluid 
inertial waves, circularly polarised. The second pair of modes 
(|37p are different, representing the proton inertial modes as 
viewed in a frame rotating with the neutrons. These modes 
are always damped by mutual friction. 

Next, we derive approximate solutions for weak and 
strong drag without any limitation on the wavelength. The 
TZ case is the easiest. Retaining only the leading or- 
der mutual friction, the modes have the form (|36p - (|37| l with 
B' = 0. As we have just discussed, these are neutron inertial 
waves weakly coupled to the protons and vice-versa. 

The strong coupling solutions are obtained after letting 
B' —> 1 and B — > in the dispersion relation. Keeping the 
same identification for the modes as in (I36I)-(I37I). we find 



0-1,2 « ±i p|i 



1 \fl 



knwii 



-!-|(nJl + ^nfxl + 4n„[(l -f 3Xp)nr, - Xpill] 

-1 1/2 

=F8a;pnnfc||Ui||| (38) 



+ ]n, 

Xp _ 



fcll Wll 



+ ^{(^1+ ^nf^l + 40.4(1 -+- 3Xp)nn ~ Xpill] 

-1 1/2 

TSxpQnkiiwA (39) 

It is clear that the strong drag coupling alters the waves in 
the system significantly. It is no longer the case that the neu- 
trons and protons can oscillate "independently". Moreover, 
the above results show that the system may now be unsta- 
ble already at leading order (i.e. when B' = 1 and B = 0). 
In order to discuss this strong drag instability further we 
assume that Qp and iT^ are almost aligned. As we will see 
later, this assumption is valid for the free precession motion. 
Setting Qp ~ Qn the mode-frequencies become 



0-1,2 



±nn (1 - — ) - k 



- — + xpf T 2n„fc||U)||a;p] '''^ (40) 



0-3,4 ^ ±Sln 



— 1 I — k\\w\\ 



+ — + Xpf =F 2n„fc||it)||a;p] '^"^ (41) 



From these results we see that the waves are unstable if 

nn{l+Xp)'^ 



(42) 



This criterion is clearly different from the weak drag result. 



psp . Since Xp is expected to be small, the critical relative 
flow needs to be larger to trigger an instability in the strong 
drag case. However, the unstable modes no longer depend on 
the mutual friction coefficient B, as they did in the weak drag 
case. Instead, the key coupling that facilitates the instability 
is due to B' ~ 1. 

It is relevant to make an observation at this point. While 
the weak drag ins tability is well known from studies o f super - 



fluid Helium, see Glaberson. Johnson fc Ostermeied (|l974h 



ISiderv. Andersson fc Comerl ( 200?! ') for discussions, we be- 



lieve that our work provides the first suggestion that an 
analogous instability may operate in the strong drag prob- 
lem. Hence, the above results are exciting. Having said that, 
we can no longer rely on the analogy with the Helium case. 
In the weak drag situation there is a natural link between 
the instability and the transition to superfluid turbulence 
(| Andersson. Siderv fc Comeij|2007l ). In the strong drag case, 
it may seem natural to expect that the unstable modes still 
lead to the formation of vortex tangles. However, given the 
stronger coupling between the two fluids in the system the 
problem is different, and will require further analysis in or- 
der to be understood. 



3.2.2 Adding shear viscosity 

The previous results were derived assuming vanishing shear 
viscosity, for the purpose of understanding the properties of 
the various types of inertial waves. However, viscosity could 
play an important role in the problem due to the small scales 
involved. At the end of the day, if an instability is triggered 
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in the inertial waves shear viscosity should be the main coun- 
teracting agent and in order to arrive at reahstic estimates 
for the growth time of the instabihty we need to account for 
it. Having said that, it is important to bear in mind that the 
various inertial waves will be affected by viscosity in drasti- 
cally different ways. The key point is that only the proton 
fluid is directly coupled to the mechanism that generates 
viscosity (recall that protons, electrons and the magnetic 
fluxtubes are all assumed to be "locked" together) while any 
communication with the neutron component is mediated by 
the mutual friction. This feature is clearly encoded in the 
hydrodynamical equations (IT|-@. 

The interplay between mutual friction and viscous cou- 
pling is exposed if we amend the approximate mode solutions 
from the preceding Section with viscosity. Considering first 
the short wavelength limit we now find 



(71,2 ~ ±20,1 + {iB T S')(2!^n T h'^w 



± 



l^cc ^ 1 1 



[B^ - [B'f] (43) 



0"3,4 



+ ^(jB^e') + ii.cofcif 



(44) 



where in each case we have kept only the leading order vis- 
cous tern|f|. The key result is that only the "proton" modes 
(73,4 are damped by viscosity in the usual way, that is, via 
fcy term which dominates at small scales. In contrast, 
the viscosity dependence of the "neutron" modes (7i,2 is in- 
timately linked to the mutual friction coupling, leading to a 
perhaps surprising scaling with negative powers of A;|| . 

As in the inviscid problem, the TZ solutions simply 
foUow from dHJ and (gH) after setting B = B' = 0. In this 
limit the (7i,2 modes are essentially not affected by viscosity. 
This is not surprising since they describe oscillations in the 
neutron fluid only. 

In the strong drag limit we find the following viscous 
modes (assuming Qp ~ Sin), which generalise (|40|l and (|4ip . 



(71,2 ~ ±fln ^1 — —J + + ^i^ccfcy 



+ nl{l + xp)'^ T'2^nki\wiiXp\''^^ (45) 



(73,4 ~ ±f2n ^1 - -I- fc||W|| -I- ^l^cckl 

+ — [ - ^{UcckflXp)'^ ±iXp{Xp - l)f2n!^Gcfcjj 



+ ^l'i{l + xpf T'i^nkuwiiXp'^''^^ (46) 



Not surprisingly these expressions admit an instability at 
leading order (top pair), but they also demonstrate the con- 
trasting viscosity dependence. To see this, consider the first 

^ It is worth noting that these results are obtained by holding 
Uee fixed and then taking the large fcy limit. In efl^ect, this means 
that one cannot use these expressions to deduce the Uee 
behaviour. 



pair of solutions and assume that the viscosity is negligi- 
ble. Then increase feefcy while keeping all other parameters 
fixed. It is then easy to see that, once the leading viscous 
term under the square-root becomes dominant, there will be 
a cancellation with the other ifeefc|/2 term. Thus, the lead- 
ing viscous damping term does not have the expected form, 
it is much weaker. This is in contrast to the result for the 
second pair of solutions. In that case the two viscosity terms 
combine to give the anticipated ivcck^^ damping term. 

AU these approximate viscous solutions convey the same 
physical information. The modes that represent pure neu- 
tron inertial waves in the limit of vanishing mutual friction 
exhibit an unconventional - and much weaker - dependence 
on shear viscosity. These are also the only modes that can 
become unstable and lead to the formation of vortex tangles 
and turbulence. 



3.3 Growth timescales 

The numerical solution of the full system (|32|l - (|35| l unveils 
the presence of an instability for a large portion of the pa- 
rameter space and, as we show below, the associated growth 
time Tgrow = 1/Im a can be quite short. The numeri- 
cal timescales can be well approximated by the preceding 
analytic mode solutions (for difi'erent wavelength ranges). 
Hence, we are able to provide simple analytic estimates for 
Tgrow . These estimates will be a useful tool in our discussion 
of precessing neutron stars. 

The first estimate follows from the short wavelength 
resuhs, (I36l or (|43|l . 



2r2nS V 



- 1 



(47) 



Based on our previous discussion we expect this result to be 
accurate also in the weak drag limit, irrespective of wave- 
length (provided of course that X <^ R). 

A second estimate is relevant in the strong coupling 
regime. Since the relevant mode solutions are sensitive to 
viscosity, we need to distinguish between the inviscid and 
the viscous case. For the former we find from (|40|l . 



Xp{l 



-1/2 



(48) 



In the viscous case it is not easy to extract a similarly simple 
estimate. This is not surprising since the modes (|45p are 
more complicated than the inviscid ones. Hence, we calculate 
1/Im (7 numerically in this case. 

Our results for the instability growth times are illus- 
trated in Figures [l][5l In the figures we compare the different 
approximations for Tgrow to the full numerical results. The 
data in the figures are obtained using a w^^ corresponding to 
a mismatch between flp and Q]^ of 1° and a rotation period 
of 1 s. Later we will relate these results to actual precession 
solutions. 

Let us first consider Figure [1] where we show the growth 
time as a function of the wavelength for fixed temperature 
and drag coefficient. The data concerns the intermediate and 
weak drag regimes, explicitly TZ = 1 and TZ — 10"''. As we 
discussed earlier, there will be a short wavelength cut-off. 
A more detailed analysis would link this to the validity of 
the averaging that led to the macroscopic hydrodynamics 
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Figure 1. Numerical results (solid lines) for the instability 
growth time Tgi-ow as a function of the wavelength A, for two 
choices for the core temperature, T = 10^ K and 10^ K (left and 
right panel, respectively). The analytic short wavelength approxi- 
mation discussed in the main text is superimposed (dashed lines). 
The results represent weak (Ti = 10"'^) and intermediate {TZ = 1) 
drag. The angle between Hn and Hp is taken to be 1° and the 
rotational period is P = 1 s. The short- wavelength cut-off, which 
is about 1 cm for this rotation rate, is indicated by a grey area. 



equations ([TJ and This is a very dilRcult problem. To 
make immediate progress, we will simply assume that the 
instability analysis becomes invalid once the wavelength is 
so short that the fluid description is no longer valid. Then 
it seems reasonable to use something like 

1/2 

A„,in ~ lOOdn ~ ( ^ 1 cm (49) 



P 
Is, 

The expected short-wavelength cut-ofT for a rotation period 
P = 1 s is then about 1 cm, as indicated in the figure. From 
these results it is clear that the short wavelength approxi- 
mation (I47|) is very good for (essentially) the entire A range 
and all relevant temperatures. 

In the strong drag regime, for TZ ^ 1, the variation of 
growth time with wavelength is more complicated, cf. Fig- 
ure[5]which shows data for TZ = 10^. The small A approxima- 
tion again matches the exact result well in the region near 
the short wavelength cut-off, but it is not accurate for all 
wavelengths. For intermediate wavelengths the strong drag 
viscous solution (|45|) is more suitable. By combining the two 
approximations we arrive at an almost perfect agreement 
with the numerical results, including the transition region 
where shear viscosity is relevant. It is also apparent from 
the figure that the inviscid estimate (|48p is accurate only 
for relatively hot stars, for which shear viscosity is weak. 
This is as expected. 

We get a different view of the results if we fix the wave- 
length and vary the drag coefficient TZ. Figure |3] shows such 
results for a rotation period of P = 1 s at the corresponding 
short wavelength cut-off, i.e. A = 1 cm. From the left panel 
we see that the small A approximation (|47|l is accurate for 
TZ < lO' or so at T = lO'^ K. The remaining parameter 
space is well approximated by the viscous 7?, — > oo solu- 
tion (|43|) . The inviscid strong-drag approximation (|48|l is 
never really adequate for this very short wavelength. Simi- 
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Figure 2. Same as Figure [T] but for strong drag Ti, = 10^. In 
this regime, shear viscosity affects the results for intermediate 
wavelengths (around A = 1000 cm and 10 cm in the left and 
right panels, respectively.) The various approximations discussed 
in the main text are superimposed (dashed lines) . It is important 
to note that the unstable modes are not strongly affected by shear 
viscosity in the short wavelength limit. 



lar results for A — 100 cm are shown in Figure U Note that 
in both figures the growth time is shortest for 71 of order 
unity, and the general behaviour follows immediately from 
Tgiow 1/B ~ 7?, + 1/TZ. This is analogous to the fact that 
the glitch relaxation timescale is exp ected to bo long in both 
the weak and the strong drag limit (lAlpar fc Sauls 1988). 

One can show that for very strong coupling, TZ > 10^ 
or so, the fastest growing waves are those with moderately 
short wavelengths, A > 100 cm. For such large values of TZ 
the mode solution (|45p is accurate even for very short wave- 
lengths. Whether this very strong drag regime is physically 
relevant is not clear. 

From these comparisons with the numerical results we 
arrive at a practical guide to the validity of the analytic 
timescale estimates. For a wide range of TZ the fastest grow- 
ing modes are in the short wavelength regime near the cut- 
off, i.e. A ~ 1 cm for P = 1 s, where the estimate (|47|l 
applies. For very strong coupling (depending on the tem- 
perature) the instability grows faster in modes with moder- 
ately short wavelengths (above A ~ 100 cm) in which case 
(|45[) applies. In fact, this approximation tends to be good 
for intermediate wavelengths at all temperatures for strong 
drag. 



4 ARE PRECESSING NEUTRON STARS 
STABLE? 

In order to address this key question we first need to discuss 
the "standard" two-fiuid model for neu tron star precession 
( Sedraki an. Wasserman fc Corded [l999l ) in some detail. A 
self-contained discussion of the relevant precession modes 
can be found in Appendix |X] We want to combine these 
precession solutions with the preceding plane-wave analysis 
to establish whether the superfiuid instability is present or 
not. In principle, this analysis may provide us with a "no-go 
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Figure 3. Numerical results (solid lines) for the instability 
growth time Tgrow as a function of mutual friction drag TZ for 
two choices of core temperature, T = 10^ K and 10® K (left and 
right panel, respectively). The analytic approximations discussed 
in the main text are superimposed (dashed lines). Here we have 
taken the wavelength to be that of the short wavelength cut-off, 
i.e. for P = 1 s we have A = 1 cm. The angle between Qn and Sip 
is assumed to be 1° 
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Figure 4. Same as Figure [S] but for A = 100 cm. 



theorem" for neutron star precession. After all, the preces- 
sion modes are obtained by assuming global solid body rota- 
tion. If the corresponding solution exhibits a rapidly growing 
short wavelength instability, then the assumption of solid 
body rotation does not hold. The precession solution would 
simply be inconsistent. The upshot of this would be that one 
cannot neglect the hydrodynamical degrees of freedom of the 
precession problem. Any analysis of the precession problem 
would have to account for the hydrodynamics of the neutron 
star core, which would be a serious computational challenge. 

In our analysis of the precession modes in 
Appendix \K\ we go beyond t he o riginal work of 
ISedrakian. Wasserman fc Corded (|l999l ) by providing 
general expressions for the precession period Ppr, the damp- 
ing time Td and the angular amplitude 9^ (the so-called 



"wobble" angle) for arbitrary mutual friction coupling. 
These are the key parameter^ which enable us to make 
contact with the two-fluid instability analysis. However, it 
turns out that the range of TZ for which precession is not 
overdamped by mutual friction (in the sense that Td > Ppi 
) is entirely cover ed by the approximate (weak an d strong 
drag) results of ISedrakian. Wasserman fc Corded ([l99^. 
This is clearly illustrated in Figure lAll Therefore we only 
discuss results for the stability of precession in the limits of 
weak and strong drag. 

From the practical point of view, the perturbative cal- 
culation in Section|3]assumes a static background. Therefore 
a given precession mode can only serve as background pro- 
vided that the unstable waves grow fast enough. To ensure 
that this is the case, we we will require Tgrow <C nnn(Ppr, Td) 
for the instability to be astrophysically relevant in a precess- 
ing star. 



4.1 Weak drag precession 

As described in Appendix|Xl there are two precession modes 
for any given value of the mutual friction parameter TZ. Let 
us first consider the case of weak drag precession {TZ ^ 1) 
and the mode that, in this weak coupling limit, reduces to 
the standard classical free precession of a biaxial body (i.e. 
to the sort of precession commonly cons i dered in the lit- 
eratur e, see e.g. IJones fc Andersson (|200ll ): lLink fc EpsteinI 
(|200ll )). This mode represents weakly damped, long-period 
precession with 



P 

e 



(50) 



where e is the (dimensionless) deformation associated with 
the biaxiality in the crust's moment of inertia (see Ap- 
pendix lAjl . For PSR B1828-11 this deformation is estimated 
to be e « 10"*. Moreover, assuming that PSR B1828-11 is 
described by the weak drag precession solution, one should 
be able to constrain B using t he fact that since its discovery 
( Stai rs. Lvne fc Shemaij[200Gl ) several undamped precession 
oscillations have been monitored. If A'' is the number of os- 
cillations then we should have Td > NPpi from which we 
find 7^ < {0.16/N){h/Is). Using N x 10 and 7i//, « 0.1 
we find TZ < 2 X 10^^. Combining this upper limit with 
the lower limit placed by Vela's glitch relaxation timescale 
we find that only a relatively narrow interval of weak drag 
mutual friction is allowed, 2 x 10~^ < TZ < 2x 10~^. Interest - 
ingly, the value suggested bylAlpar. Langer fc Saulj l|l984 ): 
lAndersson. Siderv fc Comerl l|2006h . TZ W^. falls within 
this range. 

The weak drag instability criterion (I25p implies that 
waves with wavelength shorter than 



(51) 



are unstable. For the given weak drag precession solution we 
have (Appendix |AJ 



(52) 



* Note that Ppr and 9^ are in fact observable quantities, 
extracted from the timing data of the cand idate precessors 
lljones fc AnderssonllioOll : iLink fc Epsteinll200ll) . 
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which means that 



5.4 X 10" 



1=^ 



10- 



R 



lOBcm 



(53) 



Since we need to have Amax > Amin in order to argue that 
the instability is relevant we see that we must have 

1/2 

> 1900 ( ^ 1 ( rrr-^ ] ( -r,;^ \ (54) 



VlQf'cm 



In words, this inequality gives the wobble angle required for 
the fluid in at least some portion of the star to be unstable 
for inertial waves with A = Amin. 

To demonstrate that the instability analysis is consis- 
tent we need to show that the growth time for the unstable 
modes is short compared to the precession period. Using 
(|47p . which according to Figure [1] matches the exact weak 
drag numerical results, we find 



A 



2-kB\'W\\ 



(55) 



This should be a good approximation for the shortest wave- 
length modes in the problem. Using the relevant result (|52p 
for w\\ we find that 



> 



4x10' 



3 



4 X 10- 



10- 



P_ 

Fs 



R 



On the other hand, the precession period is 

f P 



(56) 
(57) 



Requiring that rgrow < 0.1 Ppr we see that we must have 



A < 25 



4 X 10- 



1° / V 10'' cm 



(58) 



We need to ask whether there is room for unstable waves 
with wavelengths longer than Amin. This would be the case 
if 



> 4 X 10" 



B 



4 X 10-* 



R 



10*5 cm 



P_ 

Fs 



1/2 



(59) 

In words, this inequality gives the wobble angle required for 
the unstable A = Amin inertial waves in some portion of the 
star to have a growth time of less than a tenth of the free 
precession period. Compared to (|54|) this restriction on the 
wobble angle is clearly much less stringent. Hence, we learn 
that whenever the instability is active the growth time Tgrow 
is much shorter than the precession period. 

If we associate th e instability with the onset o f supe r- 
fluid turbulence (see I Andersson. Siderv fc Cornerl (120071 )) 
then only stars that satisfy the criterion (|54p are likely to 
be turbulent. Note that this hardly constrains the realistic 
parameter space since according to the data all candidate 



precessors have e ~ 10 



10' 



Hence it would seem safe 



to conclude that this instability does not play a role in sys- 
tems where the drag is weak, unless they are significantly 
deformed. However this does not mean that the instability 
criterion (|54p is not astrophysically relevant. For example, a 



millisecond pulsar with (say) P = 1.5 ms and e = 5 x 10~^ 
(which is a realistic value) would undergo an instability for 
a wobble angle 6^ > 1.5°. A more extreme example would 
be the su ggested sighting of free pr ecession of the SN1987A 
remnant l|Middleditch et ah 1 l2000al lbl). Taking a spin period 
of 2.14 ms and a free precession period of ~ 10^ s, we would 
have e ~ 2 X 10~®. For these values, the instability crite- 
rio n (I54p would be sat ifisfied for 6^ > 0.4°. As discussed 
by I Jones fc AnderssonI (|2001i ) , the wobble angle would have 
to be much larger than 1° in order to explain the observa- 
tions. In other words, this system would suffer the superfluid 
instability. 

Now consider the second mode in the weak coupling 
limit. As discussed in Appendix [X] in the case of zero cou- 
pling this mode simply corresponds to a misalignment be- 
tween the rotation axes of the neutron and proton fluids, 
with both angular velocity vectors fixed in the inertial frame, 
i.e. no precession. Weak coupling then serves to induce a 
slow precession, with the star rotating about the 3-axis of 
the biaxial crust, and with this axis itself tracing out a cone 
at the longer free precession period. This is the opposite of 
the standard free precession motion of the first mode, where 
the 3-a:xis traces out a cone in space at the "spin" frequency, 
with a slow additional precession rotation superimposed. 

However, from the observational point of view, both 
motions would produce modulations in the radio data and 
so this second mode cannot be discounted. A crucial differ- 
ence is that for this second mode the ratio Ppi/P is deter- 
mined not by the biaxiality of the crust e, but instead by 
the strength of the coupling. We have 



Ti- 



ll 



Is+Il 



(60) 



The observed P/Ppr ^ 10~* would then give an estimate 
oi TZ ^ or less. This is smaller than the typical esti- 

ma tes based on microphysical consider ations (see for exam- 
ple [Andersson^mider^^j^ome^ (|2006l )) and is also smaller 
than the lower bound placed on TZ from glitch observations 
(see discussion at the end of Section [2|, arguing against this 
mode as an explanation of the PSR B1828-11 observations. 

A further argument against this mode being relevant 
comes from our stability analysis. For this mode we find 
that 

«;|| < le^nR (61) 
and the instability is present for wavelengths shorter than 

Amax = Tvie^R (62) 

where 

h 



1 = 1 + 



(63) 



This means that, if we require Amax > Amin then an insta- 
bility is present if 



, , > 2 X 10 , ^ 
1° J I \ls 



IfPV'^f R 



106 cm 



(64) 



Alternatively, this can be expressed as a constraint on the 
rotation period. To avoid the instability we require 



P > 3 X lO^P f f — ' 

1° / V lO*" cm 



(65) 
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As one would have expected given the large relative flow 
associated with it, this mode is seriously affected by the 
superfluid instability. The growth time is 



r^row > 3.6 X lO^i ( ^ 



/ Vlsy V4 X 10 



B 



1° 



(66) 

If we require that this is smaller than a tenth of the preces- 
sion period, then we must have 



A < 3 X lO^J 



B 



(67) 



4 X 10-4 J \ l\ 

Finally, requiring this to be larger than the cut-off Amin we 
see that the instability is present for 

I > .X X 111 "^1 I 

1° 



/ V 4 X 10-* 



Is 



(68) 



The analysis suggests that this mode would be unstable for 
any signiflcant misalignment. It is therefore unlikely to be 
present in real stars. 



shown in Figure [5] The numerical timescales corroborate 
the analytic estimates. For a large part of the parameter 
space the instability grows at sub-millisecond timescales, 
much faster than the precession period. This would be the 
case, for example, for PSR B1828-11 (P f» 0.4 s, T ~ lO'^ K) 
if it were to precess under strong mutual friction. 

As we have already discussed, the strong-drag results 
are temperature dependent. This is clear from the results 
shown in Figure (5] For very large TZ, the fastest growing in- 
stability is no longer associated with the short-wavelength 
cut-off region. Instead, the fastest Tgrow is temperature de- 
pendent. This longer wavelength part of the instability win- 
dow can be understood from (|45p and (|46p . For longer wave- 
lengths (small fc||) one would expect shear viscosity to be 
unimportant. Thus the modes of the system should be well 
approximated by the inviscid result. In the inviscid case (|45|) 
implies that waves with wavelength shorter than 



A„ 



47r|TO|i 



f2n (1 + Xp) 



47ra;p|ui||| 



(74) 



4.2 Strong drag precession 

For strong mutual friction the nature of the precessional 
motion chan ges signifi c antly. To some extent, this was first 
discussed bv IShahamI (|l977l ). see Appendix 1X1 for details. 
Precession is fast, with period 



P[e+ — 



(69) 



As we discussed in Section 13.31 in the case of strong 
drag the instability is sensitive to the action of viscosity 
and as a consequence the "window" of the fastest growing 
wavelengths shifts with varying TZ. The most natural choice 
is to assume that 1 <^TZ < 10^ (see discussion at the end of 
Section[21) which means that the small A approximation (|47p 
for Tgrow describes the fastest growing waves accurately. We 
again have (I55p . which can be combined with the obtained 
approximation for (Appendix IX)! 



< 



For typical parameters we find 
W>140(|f) 



Tl_ 

103 



R 



(70) 



(71) 



For consistency Tgrow needs to be shorter than the precession 
period ([69|. Requiring Tgrow < O.lPpr as before, this restricts 
the wavelength to 



A < 70 



103 



(72) 



The short lengthscale cut-off obviously remains as in the 
weak drag case. Thus we see that we need 

1/2 

(73) 



> 1.4 X 10" 

1° 



TZ 
103 



in order for the instability to be relevant. This clearly sets 
a very severe constraint on the fast precession solution. 

The above results are, however, only relevant for short 
wavelengths. It is, of course, straightforward to find the gen- 
eral solution numerically. This provides the exact Tgrow as 



will be unstable. Working out the associated growth time 
we find 



(75) 

This estimate should be valid for long wavelengths (< Amax) 
and/or high temperatures. As we move towards shorter 
wavelengths at fixed T, the importance of shear viscosity 
will increase, cf. Figure [S] To estimate the point at which 
the viscous contribution becomes important we can balance 
the first and last terms under the square-root in (|45p . That 
is, we consider 



(76) 



This suggests that the viscosity starts to play a role when 



A 360 



Is 



10 km 



VlQSK 
(77) 

Combining this with the inviscid growth timescale (|75p we 
find 



:6x10-M'£) 



1° ) VO.l/ VlO^K ) 



(78) 

Adding the estimated (A, Tgrow) points to the data in Fig- 
ure[5]we see that we have arrived at a reasonable approxima- 
tion for the fastest growing mode in the intermediate wave- 
length regime. Finally, we can use this estimate to check 
when one would expect there to be such an instability win- 
dow. If we require Tgrow < O.lPpr as usual, then we find that 
the longer wavelength instability is present for 



T > 1.5 X l(f ( ^ 
Is 



1/2 



To conclude, our results provide a strong indication that 
the fast precession solution may not be relevant for realistic 
systems. Since it is generically unstable to small scale inertial 
waves the solid-body rotation assumption does not hold and 
the solution is inconsistent. 
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R=io' R=io'' 




10"' 10° lo' 10^ 10^' 10* 10^ 10'' lo" lo' 10^ lo' lo" 10^ 
wavelength X (cm) 

Figure 5. Instability growth times rgrow(A) for strong drag fast 
precession with wobble angle = 1°. Two choices for the drag 
coefficient are shown, Ti. = 10^ and 10^ (left and right panel, re- 
spectively). The instability is present for all wavelengths down 
to Amin (shaded area) , despite the action of viscosity. For a wide 
range of core temperatures rgrow is well below the precession pe- 
riod Ppr ~ XpP (we have used P = 1 s and Xp = 0.1) which is 
shown as a horizontal dashed line. The analytic approximations 
discussed in the main text are superimposed as dashed curves for 
the T = lO''' K case. Estimates for the minima in the T = lO'^ K 
and 10^ K curves, obtained from I I77I I and I I78II are indicated by 
-I- symbols. 



5 SUMMARY AND OUTLOOK 

Within the standard two-fluid model for superfluid neu- 
tron stars we have discussed a new class of instabilities 
that set in through short wavelength inertial waves, pro- 
vided there is sufficient relative flow along the neutron vor- 
tices. As a demonstration of the astrophysical relevance of 
these instabilities we have considered freely precessing neu- 
tron stars, in which the two angular velocity vectors are 
naturally misaligned. Our analysis is based on the usual 
form for the superfluid mutual friction, which depends on 
the strength of the drag that determines the friction be- 
tween neutron vortices and the proton/electron fluid. Our 
results show that precession may trigger instabilities both 
in the weak- and strong drag regimes. Based on the anal- 
ogy with the so-c alled Donnelly-Glaberson instability in 
superfluid Helium |Glab erson. J ohnson fc Ostermeierll 19741 : 
iBarenghi. Donnellv fc V incn 20o3), we argue that this in- 
stability is likely to lead to the formation of vortex tangles 
and generate turbulence i n superfluid neutron star interiors. 

The standard model (jSedrakian. Wasserman fc Corded 

1 19991 ) predicts two possible precession regimes. Weak drag 
leads to slowly damped long period precession, while strong 
drag implies that precession should be fast. Assessing the 
stability of each precession mode we conclude that slow pre- 
cession is generally stable. In particular, this should be true 
for the few known candidate precessors, like PSR B1828- 
11. An interesting exception could be slow precession of 
millisecond-period neutron stars. These may become unsta- 
ble above a wobble angle 9-w ~ 1.5° assuming a crust defor- 
mation e_~_5_2<.J^23j,_^£ll, within the theoretic ally allowed 
range (|Link fc Cutlerll2002l : iHaskeU et al.ll2006l ). 



In contrast, we find the fast precession solution to be un- 
stable under generic conditions. This suggests that the solid- 
body rotation assu mption tha t leads to the precession mode 
( Scdrakian. Wasserman fc Cordes.1999 ) is inconsistent. The 
upshot of this could be that fast precession is not realised in 
an astrophysical system. It should certainly be the case that 
one cannot ignore the small scale hydrodynamical degrees 
of freedom in the neutron star core in an analysis of the 
precession problem. This is an interesting conclusion since 
it adds significant complexity to the modelling. It also calls 
into question some well established i deas. 

Since the seminal discussion of lShahamI (|l977l ). it has 
been known that precession is sensitive to the neutron vor- 
tex dynamics. This is natural since the vortices provide the 
bulk of the star's rotation. If the vortices interact strongly 
with the charged component (e.g. the crust) then the gy- 
roscopic nature of the vortex array will inevitably lead to 
fast precession. This corresponding motion can either be 
described in terms of a "p inned" component affecti ng the 
usual Eulerian precession (jjones fc AnderssonI I2OO1I ) or as 
a st rong drag f ast precession solution in the framework of 
Sedrakian. Was serman fc Cordcs (1999). The observed long 
period precession of PSR B1828-11 is thus consistent only 
with weak mutual friction. 

However, one can argue that the interaction between 
neutron vortices and the (much more numerous) magnetic 
fiuxtubes in a type II pro t on superconductor ought to lea d to 
strong drag (|Saulslll989l : iRuderman. Zhu fc Chenlll998l) . If 
this leads to the core vortices being effectively pinned to the 
charged component, then fast precession should result (as in 
the Shaham model). If on the other hand, the precession mo- 
tion is able to force the vortices through the fiuxtubes then 
the motion should be highly dis sipative and there would be 
no lasting precession. This led iLinkI l|2003l ) to the sugges- 
tion that the apparent confiict with the observations can be 
avoided if the protons instead form a type I superconductor. 
There would then be no fiuxtubes and slow precession may 
be possible. However, this would not accord with the gen- 
erally held view that the protons in the outer neutron star 
core will condense into a type II superconductor. 

Our results add yet another twist to this story. We 
have demonstrated that the fast precession solution, that 
would result if the vortices and the type II fiuxtubes in- 
teract strongly, is generically unstable. Although we cannot 
claim to understand the dynamical repercussions of this, it 
seems reasonable t o assume that the fast solid-body preces - 
sion motion of the ISedrakian. Wasserman fc Corded (|l999l ) 
model will be significantly affected. If the instability leads to 
a state of superfluid turbulence, as the analogy with the He- 
lium problem suggests, then it could well be that the preces- 
sion does not survive. To explore this, a meaningful analysis 
must venture beyond the rigid-body precession model and 
account for detailed superfluid hydrodynamics. This may be 
prohibitively difficult. Anyway, it is clear that we cannot yet 
use the precession observations to draw any reliable conclu- 
sions concerning the state of the proton superconductor. 

The present work can (and needs to) be extended in a 
number of ways. Our uniform density neutron star model 
is too simplistic and one would certainly want to account 
for interior stratification etcetera. If one wants to consider 
superfiuid-superconducting mixtures seriously then the two- 
fiuid model probably needs to be upgraded into a more re- 
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alistic setup with three fluids (after relaxing the assumption 
of comoving charged particles), plus the magnetic field. The 
dynamics of such a model for neutron star cores is largely un- 
explored (Men dcU 1998), yet the additional degrees of free- 
dom could lead to new features. Our calculation was per- 
turbative. Hence, little can be said about the non-linear be- 
haviour of the unstable waves. Most likely this will be com- 
plicated, due to the small scales involved. Numerical simu- 
lations of this problem may be challenging, but should be 
encouraged. In fact they may be required if we are to under- 
stand the fate of a precessing neutron star in the strong-drag 
regime. Some of the ne c essary tools have already been devel- 
oped by [PeraitaeiaD (|2005l . |2006|). and there is analogous 
work for superfluid Helium ( Barenghi. Donnellv fc VinenI 
I2OOII ). Finally, at the level of linearised theory the present 
two-fluid analysis could be modified to consider the inner 
layers of the crust, where superfiuid neutrons coexist with 
an nuclear lattice. One may expect the nature of any super- 
fiuid instabilities in that region to be different. 

These are all interesting problems that we expect to 
return to in the future. 
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APPENDIX A: 
SOLUTIONS 



EXPLICIT PRECESSION 



We want to make quantitative contact between our local 
plane-wave analysis and the global motion of a freely pre- 
cessing two-fluid neutron star. In particular, we need to work 
out the two key quantities f2p and , which are required if 
we want to assess the relevance of the plane-wave instabili- 
ties. 

The precession of two-fluid models, accounting for 
the mutual friction coupling, has been discussed by 
ISedrakian. Wasserman fc Corded (|l999[ ). As we will discuss 
in this Appendix, it is relatively straightforward to extract 
the information we require from their analysis. They take 
as their starting point the two Euler equations HI and ([2]), 
with the mutual friction force given by ([5]). The precession 
equations are then obtained by assuming that each fluid ro- 
tates as a solid body, with the two rotation axes misaligned. 
The final equations of motion have a set of mode-solutions 
that represent free precession damped by mutual friction. 
Here we will briefiy summarise the precession solutions in 
the weak and strong drag limits, and then extract the quan- 
tities Qp and Wy that enter into the plane-wave analysis. 
Having done this, we will consider an alternative approxima- 
tion which allows us to extend our analysis to (essentially) 
the entire range of drag c oefficients. 

In the analysis of ISedrakian. Wasserman fc Corded 
()l999l ^ it is assumed that the neutron star crust is deformed 
in such a way that its moment of inertia tensor has diag- 
onal components Ji = I2 ~ 13/(1 -I- e), where e is small. 
Meanwhile, the superfiuid is assumed spherical and its mo- 
ment of inertia is Is- The calculation is then carried out in 



the crust frame, considering perturbations around a back- 
ground where Ql^ = Qp — Q'^ . Assuming that the precession 
modes have time dependence exp{p^lt) one finds that p must 
solve 



2 I 
P + 



where 



[i(l -e)+ i{B - iS')] p + e(l - B' - iB) = (Al) 



1 = 1 + 



h 



(A2) 



Once we have determined the precession period and 
the mutual friction damping timescale from the roots 
of (lAll), the eigenfunctions follow fro m equation (74) in 
ISedrakian. Wasserman fc Corded (|l999l ). This provides the 
relation 



On = 

where we have defined 

fin = fin,l + jf^n,2 



P 



p + i{l- B')+B 



sip = fip^i + ir2p,2 



(A3) 



(A4) 



and fin,! 5 foi' example, represents the ii component of the 
angular velocity of the neutrons as measured in the crust 
frame. 



Al The weak drag case 

In the weak drag case have B and B' ~ B^ . As always, 
there are two mode solutions, one which corresponds to the 
familiar Eulerian precession of a biaxial body in the limit 
B 0, and the other to a less familiar motion which has no 
analogue in one component systems. We will discuss each in 
turn, beginning with the Eulerian solution. To linear order 
in the mutual friction its frequency is given by 

^ = -"^(rfk 

and the mode amplitudes in the proton and superfiuid com- 
ponents are related by 

^"-[Ti7-(TfV0 + ^+^)]"^ 

Let us now assume that the eigenfunctions are normalized, 
i.e. choosing the precession amplitude and the initial phase, 
in such a way that 

pQt 



with a real. Then we can immediately read off 
fin,i = afie"*/^"* COs{2TTt/Ppr) 
fin, 2 = afie"*/^"* sin(27rt/Ppr) 
where the damping time is 

1 h 



Td = 



BeQ. Is 



while the precession period is 



(A7) 

(A8) 
(A9) 

(AlO) 
(All) 



Ppr = 27r/efi 
Using (|A6|) we find 

fin,i = Qfie"*/^"* [:j-^cos(27rt/Ppi.) 

- (1^7)3 + ' + 77 ) sin(2WPp.)] (A12) 
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fin. 2 = a 



S^g-t/r^ T— !— Sin(27rt/Ppr) 
LI + e 

(lIV {^^'^t) ^°«(2WJ'P.)] (A13) 



It should also be recalled that we have Q.n,z = Sip, 3 = Q.. 

We want to be able to use this solution as background 
for the plane-wave calculation. For this to make sense we 
need the rotation to be essentially constant. This is the case 
as long as we consider a timescale which is short compared 
to both the damping timescale and the precession period 
Ppr. Since B is small, the strongest constraint is provided by 
the precession timescale. Assume that t <^ Ppr we can use 



and 



fin 



1 + e 



Xl + D,X3 



(A14) 



(A15) 



One of the quantities we need for the plane-wave anal- 
ysis is Q,p, which represents the projection of the charged 
component rotation along the vortex array (along To 
work this out we need the norm of Ql. We find that 



a. 



1 



(1-^q2)1/2 

This then leads to 



1 +q2 



1 + 



1 + e 



fin 



Xl + X3 



(A16) 

(A17) 

fin in the 
plane- wave calculation. 

We also need to estimate the relative linear velocity 
along the vortex array, . We obtain this quantity from 



which shows that it is reasonable to take Q„ 



where 



Vp ni 



■ e il~Xk 



(A18) 



(A19) 



In other words, the frame independent scalar quantity that 
we need follows from 



From our mode solutions we find 



1 + e 



Xl 



Then 

(fin — fip) X h : 
and we see that 



eaO, 



(l + Q,2)l/2 



1 + 



X2 



' eaflx2 



(A20) 
(A21) 

(A22) 
(A23) 



Finally, we would like to be able to make direct contact 
with observations. This involves working out the "wobble 
angle" 9m associated with our precess ion solution. As in the 
analysis of ljones fc AnderssonI (|200ll l. we define the wobble 
angle to be the angle between the total angular momentum 
axis and the crust deformation axis. In the present case the 



latter is simply £3. Meanwhile, in the crust frame, the total 
angular momentum is given by 

J = {IinP^+IsQi)xi + {Iin^+IsQ2)x2 + {l3+Is)nx3 (A24) 

Taking the modulus of this, we then find the wobble angle 
from 



J ■ X3 — {I3 + Is)fl = J cos 6u 
In the weak drag case we then have 



J ~ (Is + /3)fi 



1 



2(l + e) = 



Hence 



1 + e 



(A25) 



(A26) 



(A27) 



This completes our analysis of the motion associated with 
(IMl) . 

We now turn to the second mode found for weak mutual 
friction. Its frequency is given by 



p = —i — B 



I + e 



(A28) 



while the superfiuid and proton mode amplitudes are related 
by 



fin 



iS^fin 



(A29) 



To gain insight as to what sort of motion this mode corre- 
sponds to, consider the case B — 0. Then we can write 



fip 

fin 



0, 

afi, 



(A30) 
(A31) 



i.e. the crust is unperturbed, and keeps spinning at a rate 
fi about the inertial 3-axis, but the superfiuid does some- 
thing more complicated. Writing the superfiuid components 
explicitly, we have 



fin,i = Qficos(fit), 
fin, 2 = — Qfisin(fit). 



(A32) 
(A33) 



It is then straightforward to show that the rate of change of 
fin with respect to the inerttal frame is zero: 



dfi„ 



dfi 



dti 



dt 



"1" ^ijk^ ^ ^ 'n 



0, 



(A34) 



where the subscripts T' and 'R' label time derivatives in the 
inertial and rotating frames, respectively. It follows that the 
superfluid's angular velocity is fixed in space; in this limit 
the mode corresponds to a misalignment between the crust 
and superfiuid spin axes by an angle a. The relative velocity 
between the crust and superfiuid is then 



^ijk^nX , 



(A35) 



which gives a relative flow along the vortices at time t = 
of 



■ X2CtO,. 



(A36) 



This remains true to leading order in the case of small but 
non-zero B. 
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fin, 2 
fin 



'p,l 



fi 



P,2 



When S 7^ we instead have 
aficos(fi^)e"*''^^ 
-afisin(fi^)e"''^^^ 

QfiB^ sm(fif)e~*/"'' , 
h 



(A37) 
(A38) 

(A39) 
(A40) 

Bfi/>.- 

The total angular momentum is the sum of the neutron and 
proton contributions, which can be shown to be given by 

B 



where the damping time is given by 

1 1 + e 

Td 



Jl = 



Qffie 



/. cos 



nt - 



— Qffie '''^''JaSin 
/sfi. 



fit - 



1 + 



l + £ 



(A42) 

(A43) 
(A44) 



From these equations we can easily compute the angle be- 
tween Ji and the crust's 3-axis: 

cos-^(J-f3) = Qe"*/"" r ^\ ■ (A45) 
Is + h 

The constancy of this angle (aside from its slow monotonic 
decrease) shows that the crust's 3-axis moves on a cone 
about the fixed total angular momentum. To see exactly 
how the body moves it is best to make use of Euler angles 
{9, (f), tp) giving the orientation of the body axes with respect 
to some fixed inertial axes. Then the angular velocity com- 
ponents of the cr ust referred to the rotating crust frame take 
the form, see e.g. iLandau fc Lifshitd l| 19761 ). 

fip,i — (j) sin 6 sin tjj + 9 cos tl) , (A46) 
fip,2 ~ (j) sin 6 cos tp ~ 9 sinil) , (A47) 
fip,3 = <}}cos9 + ip. (A48) 

In defining {9, (j), ip) we are free to make any choice of fixed 
inertial frame we want, but obviously it would be simplest to 
choose one where the inertial 3-axis of this system lies along 
the fixed angular momentum vector. Then 9 — cos~^ ( J ■ X3) 
which is a constant, so the above equations reduce to 



fip,i — (j) sin 9 sin tp , 

fip,2 ~ (j) sin 9 cos tp , 

fip,3 = (j) cos 9 + 1p, 
where 



Now compare with 
fip.i — A sin fit, 
fip,2 — yl cos fit, 
fip,3 ~ fi, 
where 



+/3 



A = afie 



-t/^di 



I3' 



It follows at once that 



= fiS 



Is+h 

h 



(A49) 
(A50) 
(A51) 

(A52) 

(A53) 
(A54) 
(A55) 

(A56) 
(A57) 



and the fip.s equation gives 



V> = fi - 



fi 



1-B 



Is+h 



h 



Collecting results: 



-t/^d 



Sfi 



Is+h 
Is+h 



h 



fi 



l-B- 



h 



h 



constant, 

< fi, 
^ fi. 



(A58) 

(A59) 
(A60) 
(A61) 



We therefore see that the motion is still rather simple, with 
the crust's 3-axis moving in a cone of half-angle 9 about 
the fixed angular momentum axis at the slow rate tp <^ Q, 
with a rapid rotation at rate ~ fi about the crust's 3-axis 
superimposed. 

Amusingly, this is the opposite of normal free preces- 
sion, i.e. small angle free precession of a single component 
rigid biaxial (but nearly spherical) body, as in that case 
« fi and li/i] « efi ^ fi. However, from the point of view 
of the radio observations, both sorts of precession consist 
of two superimposed rotations misaligned by a small an- 
gle and so are observationally indistinguishable. This means 
that rather than interpreting PSR1828-11 as undergoing the 
usual precession mode with the ratio P/Ppr ~ e, we would 
instead have 



^^B 

rP 



Is+h 

h 



(A62) 



The potential physical significance of this is discussed in 
Section |4T] 



A2 The strong drag case 

In the strong drag limit we have B 1 and 1 — S' ~ . In 
this case we find the two precession modes 



^[i{l-B')+B] 



(A63) 



a mode which damps out before it completes one full cycle, 
and 



P ■ 



B 



1 + aIs 



In these expressions we have used 



(A64) 



(A65) 



The solution (|A64|I represents slowly damped fast preces- 
sion, since h/Ii '-^ 10 ^ e. This is essentially the precession 
mode discovered byHhaha^ Il973) for a neutron star model 
where an amount h of neutron vortices is perfectly pinned. 
Focusing on (|A64|) we find that the damping time is 



Td 



h 



(A66) 



Sfi(l + a) Is 
while the precession period is 

P 

' an 

As in the weak coupling case the damping is relatively slow. 



(A67) 
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The amplitudes are related by 

fin = -i-t^p = -iB (e + y] (A68) 



and 



Again assuming that the plane-wave analysis is relevant on 
a timescale short compared to Ppr, we find the approximate 
mode solution 



fin 



In this case we see that 



I, 



h 



h 

and it immediately follows that 

fip ~ fin 



X2 + fiis 



X2 + X3 



(A69) 
(A70) 

(A71) 
(A72) 



as before. 

To work out the relative flow along the vortex array, we 
take the same steps as in the weak drag case. This leads to 



fin — fip ~ — afiii — olBQ, ( e + 

and 

(fin — fip) X n ; 



which means that 



-aBVt ( e + 



Xl 



X2 



(A73) 



+ afii2 - a Sfi ( e + 7^ 



x-i (A74) 



(A75) 



The final step again concerns the wobble angle. For the 
above solution we find that 



J ~ (/. + /3)fi 

which leads to 



1 + 



„2 r2 
a ii 



2(/s + hY 



all 
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(A76) 



(A77) 



since e <C Is/Ii 



A3 The "general" case 

So far, the precession solutions that we have dis- 
cussed are identical to the lim it ing c ases considered by 
ISedrakian. Wasserman fc Corded l|l999l ). However, we can 
do better than this and find solutions that are valid across 
the range of permissible drag coefficients, TZ. The starting 
point is (|A1|I as before, but now we assume that e is the 
small parameter. This approximation is generally valid as 
long as B >> e, i.e. it fails only in the extreme weak drag 
limit (which we have already analyzed above). 

Considering (|A1|I to order e we find the two solutions 



Pi 



1-B' -iB 
^i{l-iB') + iB 



(A78) 



P2 



-IB -ill -IB' 



iB + B' 



(A79) 

hi{l- IB')+IB 

It is quite straightforward to show that these solution behave 
in the anticipated way in the weak- and strong drag limits. In 
particular, we find that HA78|) limits to HA5|) when TZ 
and to (|A63|) for 7?, oo. Meanwhile, (|A79P reduces to 
(|A28|) and to HA64|) in the respective limits. 

To extract the precession period and the damping time 
we need the real and imaginary parts of p. The first solution 
(|A78I) can be written 



Pi 



|-(1 -i)B-i [JB^ + (1 - S')(l 



IB 



')]} 

(A80) 



where we have introduced 

V = Pb"" + IB'f (A81) 

Since we assume that the perturbations behave as ef"' we 
see that the damping time is 



V 



Td •' 



(1 - l)e5fi 



while the precession period is 
27r 



V 



(A82) 



(A83) 



efi/S2-h(l-S')(l-/S') 

Given the above estimates we can determine the pa- 
rameter range for which this precession mode is relevant. 
A reasonable criterion would be that the precession motion 
does not damp out before a single cycle is completed. This 
means that we must have r^/Ppr < 1. For the solution given 
above we find that this is true provided that 



0.016 



(A84) 



for our canonical parameters, cf. the left panel of Figure IXT] 
For larger values of the drag coefficient this solution is over- 
damped. Since the maximum value of IZ that we need to 
consider is small, one would expect the weak drag estimate 
to be sufficiently accurate for our analysis. As we can see in 
Figure \K\\ the weak drag solution matches the exact result 
well even before TL enters the range (|A84|) . 

Let us now move on to the second precession solution, 
(|A79|) . Here the algebra gets a little bit messier. The fre- 
quency is given by 

P2^-IB+^^-^-i[{\-iB;)-e^\lB' -(1 



from which we see that the precession period is 



IB')B'^ 
(A85) 



Pp. 



27r 



1 + 



e/, IB' - (1 - IB')B' 



(1 -/B')fi 
Meanwhile, the damping timescale is 

1 



\-IB' 



1 + e- 



(A86) 



(A87) 



JfiB V hlT), 

Again, we can assess the relevance of this mode by 
checking that it is not overdamped. The behaviour is illus- 
trated in the right panel of Figure lAll In this case we find 
acceptable solution for 



7^ > 27r 



7 



(A88) 
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Figure Al. The solutions {PpnTd} for the precession modes pi 
(left panel) and p2 (right panel) for a general mutual friction drag 
7?.. These are compared to the weak and strong drag solutions, 
left and right panel respectively, represented by the dashed curves. 
The shaded areas correspond to the intervals of "R. for which pre- 
cession is not overdamped by mutual friction, i.e., r^j > Ppr. Note 
that all timescales are normalised with the rotational period P. 



In this case TZ is sufficiently large (compared to unity) that 
one would expect the strong drag approximation to give re- 
liable results. Indeed, it is clear from Figure lAll that one 
would expect the strong drag approximation to be valid 
when this mode is oscillatory. In addition, the figure shows 
that the mode is also relevant for IZ < 0.01. In that regime 
it is well approximated by the weak-drag solution HA28|) . 

In practice, the results for the general precession modes, 
Pi and P2, shows that the (much simpler) approximate weak 
and strong drag solutions apply for the entire range of TZ for 
which precession is not overdamped. 
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